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Abstract. We present preliminary results of the STAR experiment at RHIC 
on exotic particle searches in minimum bias Au+Au collisions at ^Jsnn = 200 
GeV. We observe a narrow peak at 1734 ± 0.5 ± 5 MeV in the KK'^ invariant 
mass with width consistent with the experimental resolution of about 6 MeV 
within the errors. The statistical significance can be quantified between 3 and 
6 a depending on cuts and methods. If this peak corresponds to a real particle 
state it would be a candidate for the iV" or the 'EP 1—1/2 pentaquark states. 
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1. Introduction 

Hadrons made by 4 quarks and one antiquark, called pentaquarks, were predicted 
long time ago (e.g. [T). Pentaquark states have been searched for in the past 
without success until the recent finding of a candidate for the &~^ state first by the 
LEPS collaboration [|2| and by other experiments [|S|, motivated by [0]. Candidates 
for other pentaquarks have been presented recently, in particular for the 2 (1862), 
S-(1850), S°(1864) [~ and the e°(3099) ["6". In this article we present preliminary 
results of the STAR (Solenoidal Tracker At RHIC) experiment at RHIC (Relativistic 
Heavy Ion Collider) on a search for the S" 1=1/2 as well as for the iV" pentaquark 
states in the decay mode AK^. Since there may be many pentaquark multiplets 
one may expect more than one S*^ and N'^ state to appear. 

2. Experimental Setup 

This analysis presents results from data of minimum bias Au+Au collisions at 
y^sjvjv=200 GeV, which were recorded in the 2001 run with the STAR detector 
at RHIC. A detailed description of the STAR experimental setup can be found in 
reference [E]. The present analysis is based on charged particle trajectories mea- 
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sured and identified with the help of a large cylindrical time projection chamber 
(TPC) ["S" with full azimuthal coverage, located inside a 0.5 Tesla solenoidal magnet 
allowing for momentum reconstruction. The TPC allows for the direct identifica- 
tion of charged particles with small momenta by measuring their ionization energy 
loss (dE/dx) ["S". Indirect identification of charged or neutral particles decaying at 
least partly into charged measured particles inside the TPC can be obtained with 
methods based on the topology of their decay (e.g. A, S, fl, K^ etc decays) [E]. 

3. Data Analysis Techniques 

In the present analysis we investigate the invariant mass of Ai^^ . We required the 
Z position (Z is along the beam direction) of the main interaction vertex to be ± 25 
cm around the center of the TPC. The number of Au+Au events after this cut is 
1.65 10^ and comprise the total available statistics of minimum bias Au+Au data 
at 200 GeV prior to the 2004 STAR run. We search for T^O topologies which are 
candidates for the decays A — > pTr~ (A -^ p7r+) and K^ -^ tt+tt" requiring the 
finding of a secondary VO vertex, at least 6 cm away from the primary interaction 
vertex. The Distance of Closest Approach (DCA) between the positive and the 
negative VO tracks (so called VO daughters) is required to be less than 0.8 cm. 
Each of the daughters is required to not originate from the primary interaction 
vertex (PV). In particular the DCA of the positive track to PV must be greater 
than 1 cm for the A hypothesis, 2.5 cm for the A hypothesis and 1.3 cm for the 
K'^ hypothesis. The DCA of the negative track to PV must be greater than 2.5 cm 
for the A hypothesis, 1 cm for the A hypothesis and 1.3 cm for the K'^ hypothesis. 
Since the pentaquarks are expected to decay strongly we require the yO's to point 
back to the primary vertex, in particular to have a DCA to PV less than 0.4 cm. 
We require at least 15 hits (out of maximally 45) for each track. We require that 
the A and K^ candidates which are combined to estimate their invariant mass do 
not share any tracks. The dE/dx of the VO daughters is required to be within 3 a 
around the expected dE/dx value for the assumed particle hypothesis. In order to 
enable the dE/dx identification we restrict the momenta of each track in the region 
in which a good dE/dx identification is possible [|H|, namely below 0.7 GeV for 
(anti)protons and below 0.5 GeV for pions. The latter cut on the momenta of the w 
coming from the decay of the K^ is motivated independently by Monte Carlo results 
for particles with mass in the range below 2 GeV and above the threshold for decay 
to AK'^, which show that the momentum of the K'^ is contained below 1-1.5 GeV, 
depending on the assumed initial transverse momenta of the parent. We select the 
Kg and A candidates within a window of ± 35 MeV and ± 10 MeV around their 
mean invariant masses which are shown in figure ^ (left and middle) . We accept 
only unambiguouisly identified K^, A and A as illustrated in figure^ (right) for K^ 
and A. 

4. Results 

Figure [5] shows the invariant mass of A and K^ which have been preselected as 
discussed in the previous section. The data are minimum bias Au+Au collisions 
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Fig. 1. Figures left and middle: Invariant masses of m(7r"'"7r~) and m(j)TT~) for the 
selected K'^ and A candidates. Figure right: Armenteros plot for the selected A and 

K^ candidates. 

at y/SNN =200 GeV while the upper r^ 10% of the crtot,Au+Au has been excluded 
from the analysis to reduce the background from the highest multiplicity events. 
The line shows the mixed event background expectation, which has been calculated 
using A and K^ originating from different events. 
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Fig. 2. Invariant mass distribution A K^ (in GeV) in min. bias Au+Au collisions at 
y/SNN — 200 GeV measured with the STAR experiment together with the estimated 
background using the mixed events technique (line). We use 3 MeV bins. The upper 
~ 10% of the atot, Au+Au has been excluded. 

The mixed event distribution has been normalized to the 'signal+background' dis- 
tribution in the region below 1.7 GeV and between 1.76 and 2 GeV. Other choices 
of smaller normalization regions did not influence the results significantly. We use 
bins of 3 MeV. We observe a peak at 1734 MeV. When we fit this peak with a 
Gaussian distribution restricted to the region of ± 3 MeV around the mean, plus a 
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Fig. 3. Invariant mass distribution A K'^ (in GeV) in min. bias Au-t-Au collisions at 
y^SNN = 200 GeV measured with the STAR experiment together with the estimated 
background using the mixed events technique (line). The upper ^ 10% of the 
cTtot, Au+Au has been excluded. We use 1.5 MeV bins. 



polynomial function for the background from 1.65 to 1.8 MeV, using 1 MeV bins, 
we obtain a gaussian width of 4.6 ± 2.4 MeV and a x^/DOF = 1.09. 
We obtain similar results when using a Breit Wigner distribution and different bin 
sizes. The width is consistent with the experimental resolution within the errors. 
The latter has been estimated with Monte Carlo generated particles with mass 1730 
MeV, flat dn/dy distribution ±1.5 units around midrapidity, exponential spectrum 
in rriT with inverse slope 250 MeV and Breit Wigner width of 1 MeV, decaying into 
AK^, which have been tracked through the Geant STAR simulation and have been 
embedded in real p+p STAR data. We could reconstruct the initial Monte Carlo 
particles with a mass of 1729 ± 0.7 MeV and a width of 6.3 ±1.7 MeV which is 
dominated by the experimental resolution. 

The gaussian fit to the data described above results in a mass of 1733.6 ± 0.5 MeV 
with a systematic error of ~ 5 MeV. The latter has been deduced from the deviation 
of the K'^ mass from the mean K'^ mass given by the PDG [^j in the appropriate 
region of transverse momentum. The values for the mean position of the peak and 
the width obtained by using Breit Wigner fits and/or Gauss fits and/or diiferent 
binning and/or different analysis cuts were consistent with the above values within 
the statistical errors. Assuming the mass of 1734 MeV and considering the region of 
± 1.5 a around the mean of the mass we obtain a significance of S/^/B = 30.6/V35.4 
= 5.1, while a more conservative estimate gives S/^/S + B = 3.8. The ratio of the 
signal to its own error is S/a{S) = S/\/S + 2B = 3.04. Furthermore we observe a 
peak at 1693 ± 0.5 (stat) ± 5 (syst) MeV with a significance of S/VB = 2.92. The 
width is similar to the previously discussed peak at 1734 MeV. 
Fig. |3] shows the same data as fig. |2 however with bin size 1.5 MeV, in order to 
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Fig. 4. Invariant mass distribution A K^ (in GeV) in min. bias Au+Au coUisions 
at ^/SNN = 200 GeV measured with the STAR experiment (red) together with the 
estimated background using the mixed events technique (hue). We use 3 MeV bins. 

illustrate that the 1734 peak in fig. [21 is not assymetric, but the left shoulder seen 
there is a well separated peak at mass 1726.6 ± 0.5 MeV. 

Figure 21 shows the same distribution as figure (21 but without any restriction on the 
event centrality. The significance is Sj^fB — 40.55/\/83.45 = 4.44 in this case. We 
obtain the best significance of S/^/B = 19.4/\/10.6= 5.93 for semiperipheral events. 

In the following we discuss some possible sources of systematic errors. One source of 
systematic errors especially in high multiplicity events is the fact that the tracking 
software may split a real trajectory into two tracks which will have a very similar 
momentum and can lead to selfcorrelations and to peaks in invariant mass distri- 
butions. We excluded split tracks from the analysis by demanding at least 25 hits 
on each track, as the maximum number of hits is 45. We also studied the effect on 
the invariant mass Ai^° if one selects either the same track taken e.g. as pion of 
the A and of the K^ or if one selects only the (same charge) tracks with momentum 
differences below (respectively above) e.g. 100 MeV. The 100 MeV is well above 
the expected momentum difference on the basis of the momentum resolution in this 
momentum range [|H|. None of these studies gave a peak near 1734 or 1693 MeV. 
Another possible source of systematic errors is the possibility that the proton of the 
A and the tt" of the AT^ or the A come from a A(1232) decay. As the sum of the 
A(1232) mass and the A:^(497.7) mass gives 1729.7 MeV this could lead to a peak 
at 1734 MeV, given the STAR systematic error of about 5 MeV. This possibility has 
been investigated either by cuting harder on the DCA of all tracks to the PV (as 
the A decays right at the primary vertex) or by cuting out cases giving an invariant 
mass for p{K) -\- it~ [K^) outside (respectively inside) a window of ± 30 MeV around 
the A(1232) mass. These studies showed that the peak at 1734 MeV cannot be 
understood as due to the A(1232)-|-A'g(497.7) mass reflection. 
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5. Discussion 

In this section we discuss the possible assignment of the observed peaks, if they 
would correspond to real particle states. The already known particles decaying into 
AK^ with mass near 1734 MeV namely the N(1710) and N(1720) (with Breit Wigner 
masses in the range up to 1740 MeV) have a width of 100 MeV or greater and are 
therefore not good candidates for the narrow peak seen here at 1734 MeV. This 
peak is not a candidate for the decay of the 'EP pentaquark with isospin 3/2, as the 
latter does not decay into AK° due to isospin violation. It is a possible candidate 
for two pentaquark states: the A'^'' with quark content udsds decaying into AK^ 
and the S° with isospin 1=1/2 and quark content udssd decaying into AA'''. The 
N^ can be from the antidecuplet, from an octet [El or an 27-plet [E], while the 
'EP 1=1/2 from an octet. In [[HI the best estimate of m{N) is ~ 1730 MeV. The 
mass of the N° is expected in the approx. range 1650-1780 MeV [[El CH Hg. 
The E!° 1=1/2 is expected to be near 1700 [^^ o^ i* maybe degenerate with the 
S°(1860) 1=3/2 [fnj. The fact that we don't observe a peak above background near 
1850 or 1860 MeV disfavours the latter possibility. However the branching ratio of 
a possible S''(1850— 1860) 1=1/2 to AAT^ may be small and we need more statistics 
to observe it. The mass of the peak at 1734 MeV is in very good agreement with 
the A^ mass of ^ 1730 MeV suggested by Arndt et al [E|- In this paper a modified 
Partial Wave analysis allows to search for narrow states and presents two candidate 
N masses, 1680 and/or 1730 MeV with width below 30 MeV. 

We don't yet observe a pronounced peak at 1734 or 1693 MeV in the anti-channel 
AK^. This is work in progress. If pentaquarks would be primarily formed through 
quark coalescence rather than through hadronic interactions in Au-|-Au collisions 
at RHIC, the expected antipentaquark to pentaquark ratio is estimated as follows: 
EP/EP - ^M - I . (1)2 o. 0.90 respectively A^/Af° ~ g§ ~ (^f - 0.73 while 
we used a ratio of p/p = 0.73 [ 15 and assumed J/ s =1. These values are at 
production, while further reduction of the ratios can follow from e.g. absorbtion of 
the decay products. Therefore, the non-observation of an antiparticle may favour 
the N^ hypothesis. The possible peak we observe at 1693 ± 0.5 MeV is a candidate 
for the state S(1690) with mass 1690 ± 10 MeV and width below 30 MeV [ 10,. If 
so, we improve the PDG limit of the width from F < 30 MeV at present, to F < 6 
MeV. 

6. Conclusions 

We present preliminary results of the STAR experiment at RHIC on a search for 
the A^*^ and S° 1=1/2 pentaquarks in minimum bias Au+Au collisions at 200 GeV 
through the decay channel AK'^. We observe a peak at a mass of 1733.6 ± 0.5 (stat) 
± 5 MeV (syst) and width consistent with the experimental resolution of about 6 
MeV within the errors, which gains in significance when restricting the upper ~ 
10% of (Jtot- We obtain an estimate of the significance of S/\/B = 30.6/-\/35.4 
= 5.1, S/yJWTB = 3.8 and S/a{S) = 3.04. The best significance of S/^/B = 
19.4/v'10.6= 5.93 is obtained for semiperipheral events. Systematic studies suggest 
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that this peak is not due to missidentifications or spht tracks or to the decay of 
A(1232). If this peak corresponds to a real particle state it would be a candidate for 
the 7V° (octet, antidecuplet or 27-plet) and S° 1=1/2 (octet) pentaquark states. We 
don't observe a peak near 1850 - 1860 MeV, disfavouring the picture of degenerate 
octet and antidecuplet even though a low branching ratio to KK^ may prevent us 
from observing a peak. The non observation of a possible antiparticle, while this is 
work in progress, may favour the N'^ hypothesis and dominant production through 
quark coalescence. The ambiguity between 'EP and N'^ can be resolved by searching 
for their isospin partners iV"*" -^ KK^ and S~ — > KK~ . This is work in progress 
as well as searches for the 0+ [El- Furthermore, new data taken in our 2004 run 
will enhance the statistics of minimum bias Au+Au events by a factor ^ 10-15. 
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